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Abstract. Ultracold neutral plasmas are formed by photoionizing laser-cooled 
atoms near the ionization threshold. Through the application of atomic physics 
techniques and diagnostics, these experiments stretch the boundaries of traditional 
neutral plasma physics. The electron temperature in these plasmas ranges from 
1-1000 K and the ion temperature is around IK. The density can approach 
10 11 cm -3 . Fundamental interest stems from the possibility of creating strongly- 
coupled plasmas, but recombination, collective modes, and thermalization in 
these systems have also been studied. Optical absorption images of a strontium 
plasma, using the Sr+ 2 Si/2 — ► 2 Pi/2 transition at 422 nm, depict the density 
profile of the plasma, and probe kinetics on a 50 ns time-scale. The Doppler- 
broadened ion absorption spectrum measures the ion velocity distribution, which 
gives an accurate measure of the ion dynamics in the first microsecond after 
photoionization. 



1. Introduction 

The study of ionized gases in neutral plasmas traditionally spans temperatures ranging 
from 10 16 K in the magnetosphere of a pulsar to 300 K in the earth's ionosphere. 
Ultracold neutral plasmas pQ, formed by photoionizing laser-cooled atoms near the 
ionization threshold, stretch the boundaries of these studies. The electron temperature 
ranges from 1-1000K and the ion temperature is around 1 K. The density can approach 
10 11 cm~ 3 . 

Fundamental interest in these systems stems from a range of phenomena in the 
ultracold regime. It is possible to form strongly-coupled plasmas |2j - systems in 
which the electrical interaction energy between the charged particles exceeds the 
average kinetic energy. This reverses the traditional energy hierarchy that underlies 
our normal understanding of plasmas based on concepts such as Debye screening 
and hydrodynamics. Strongly-coupled plasmas appear in exotic environments, 
such as dense astrophysical systems 3 , matter irradiated with intense-laser fields 
01 E|, dusty plasmas of highly charged macroscopic particles [2], or non-neutral 
trapped ion plasmas that are laser-cooled until they freeze into a Wigner crystal. 
Ultracold plasmas are excellent systems in which to study strong-coupling achieved at 
relatively low-density in neutral systems where recombination, collective modes, and 
thermalization are all interesting things to explore. 

Early experiments with ultracold neutral plasmas used charged particle detection 
techniques to study methods and conditions for forming the plasma .1], excitation and 
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Figure 1. Strontium atomic and ionic 
energy levels involved in the experiment, 
with decay rates. (A) Neutral atoms are 
laser cooled and trapped in a magneto- 
optical trap (MOT) operating on the 
1 5o — 1 Pi transition at 460.9 nm, as 
described in \YA\ . Atoms excited to the 
1 Pi level by the MOT lasers are ionized 
by photons from a laser at ~ 412 nm. 
(B) Ions are imaged using the 2 Si 12 — 
2 Pi/2 transition at 421.7 nm. 2 Pi/2 
ions decay to the 2 £>3/2 state 7% of the 
time, after which they cease to interact 
with the imaging beam. This does not 
complicate the experiment because ions 
typically scatter fewer than one photon 
during the time the imaging beam is on. 

detection of plasma oscillations [8], dynamics of the plasma expansion 8 , and collisional 
recombination into Rydberg atomic statesj^. In related experiments, researchers at 
the University of Virginia and Laboratoire Aime Cotton, Orsay, France studied the 
spontaneous evolution of a dense, cold cloud of Rydberg atoms into a plasma |TO] . 

Optical absorption imaging and spectroscopy, demonstrated using the Sr + 
2 Si/2 — * 2 P\/2 transition in a strontium plasma, opens many new possibilities. 
Images depict the density profile of the plasma, and the Doppler-broadened absorption 
spectrum measures the ion velocity distribution. Both can probe ion dynamics with 
50 ns resolution. This paper will concentrate on the physics we have learned from 
plasma imaging. A review of the results from studies using charged particle diagnostics 
can be found in |12| . 

Section [5] provides an overview of the creation of an ultracold neutral plasma 
and the absorption imaging technique. Section discusses the dynamics of the ions 
during the first few microseconds after photoionization. Section 01 describes the ion 
absorption spectrum, and Sec. describes recent studies of plasma equilibration. 

2. Experimental Overview 

The recipe for an ultracold neutral plasma starts with laser-cooled and trapped neutral 
atoms 14 . Photons from properly arranged laser beams scatter off the atoms to 
generate the forces for cooling and trapping. Alkali atoms, alkaline-earth atoms, and 
metastable noble gas atoms are the most common atoms for these experiments because 
they posses electric-dipole allowed transitions at convenient laser wavelengths. The 
normal trapping configuration is called a magneto-optical trap (MOT) (Fig. and |2J 

In the experiments described here, the neutral atom cloud is characterized by a 
temperature of a few mK and a density distribution given by n(r) = noexp(— r 2 /2a 2 ), 
with a ~ 0.6 mm and no ~ 6 x 10 10 cm~ 3 . The number of trapped atoms is typically 
2 x 10 8 . These parameters can be adjusted. In particular, turning off the trap and 
allowing the cloud to expand yields larger samples with lower densities. 

To form the plasma, the MOT magnets are turned off and atoms are ionized with 
photons from the cooling laser and from a 10 ns pulsed dye laser whose wavelength 
is tuned just above the ionization continuum (Fig. [IJ. Up to 30% of the neutral 
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Figure 2. From 11 . Ex- 
perimental schematic for stron- 
tium plasma experiments. The 
MOT for neutral atoms consists 
of a pair of anti-Helmholtz mag- 
netic coils and 6 laser-cooling 
beams. Atoms from a Zeeman- 
slowed atomic beam enter the 
MOT region and are trapped. 
1 Pi atoms are ionized by the 
photoionizing laser. The imag- 
ing beam passes through the 
plasma and falls on a CCD 




Figure 3. Optical depth of an 
ultracold neutral plasma. The 
delay between the formation of the 
plasma and image exposure is 85 ns. 
The plasma contains 7 X 10 7 ions 
and the initial peak ion density is 
Hoi = 2 X 10 10 cm -3 . Resolution 
is about 65 fim, limited by pixel 
averaging to improve the signal-to- 
noise ratio. 
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atoms are ionized, producing plasmas with a peak electron and ion density as high as 
noe ~ n>Qi ~ 2 x 10 10 cm~ 3 . The density profiles, n e (r) ss rii(r), follow the Gaussian 
shape of the neutral atom cloud. 

Because of the small electron-ion mass ratio, the electrons have an initial kinetic 
energy approximately equal to the difference between the photon energy and the 
ionization potential, typically between 1 and 1000 K. The initial kinetic energy for 
the ions is close to the kinetic energy of neutral atoms in the MOT. As we will discuss 
below, the resulting non-equilibrium plasma evolves rapidly. 

To record an absorption image of the plasma, a collimated laser beam, tuned 
near resonance with the principle transition in the ions, illuminates the plasma and 
falls on an image intensified CCD camera. The ions scatter photons out of the laser 
beam and create a shadow that is recorded by an intensified CCD camera. The optical 
depth (OD) is defined in terms of the image intensity without (Ibackground) and with 
(Ipiasma) the plasma present, 

OD(x,y) = \n(I backg 

round 

(x, y)/I p i asma 

(*,*/))■ (i) 

Figure 01 shows a typical absorption image. Section 0] describes how detailed 
information about the plasma is extracted from the optical depth. 
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3. Ion Dynamics 

In order to understand the details of the image analysis, it is necessary to understand 
the dynamics of the plasma. The imaging probe is most sensitive to the ion 
dynamics, so we will concentrate on this topic. The behavior of electrons was studied 
experimentally in ^ |SJ |5] and theoretically in ^1 ^5] E| • 

Ions are created with very little kinetic energy, but their initial spatially 
uncorrelated state possesses significant Coulomb potential energy compared to the 
regular lattice that represents the ground state of the system ^3[7]. As ions equilibrate 
and correlations develop, the kinetic energy increases. This process is called disorder- 
induced heating, or correlation heating, and it has been discussed in many theoretical 
papers. Early interest was generated by non-equilibrium plasmas created by fast-pulse 
laser irradiation of solid targets, for example ^3 E3 IS E2| , and the problem has been 
revisited in the context of ultracold neutral plasmas |15l llfil 1^31 124| . 

Qualitatively one expects the ion temperature after equilibration to be on the 
order of the Coulomb interaction energy between neighboring ions. A quantitative 
analysis |23| , assuming complete initial disorder and incorporating the screening effects 
of the electrons, predicts an ion temperature of 

Here, k — a/Xn where Ad = (eofesTg/nee 2 ) 1 / 2 is the Debye length. The quantity 
U = jv-e^/Lrena * s ^ ne excess potential energy per particle in units of e 2 /47reoa, 
where a = (A-rmi/S) _1 / 3 is the Wigner-Seitz radius, or interparticle distance. N. t 
is the number of ions. U has been studied with molecular dynamics simulations 
.25, for a homogeneous system of particles interacting through a Yukawa potential, 
4>{r) — 47 ^ Eor exp(— r/Aj>), which describes ions in the background of weakly coupled 
electrons!- 

For typical strontium plasmas discussed here, n w 0.1 — 1, and Arj w 2 — 8/Ltm. 
U ranges from —0.6 to —0.8, so T; is close to e 2 / 4:ire akB as expected, k is related to 
the Coulomb coupling parameter for electrons, r e , through k — \/3r e . A system 
of particles with charge q is strongly coupled when T = q 2 / AneoakBT > 1 
r e w 0.1 — 0.5 for the systems studied here, so the electrons are not strongly coupled. 
This avoids excessive complications that arise when T e approaches or initially exceeds 
unity, such as screening of the ion interaction JH] , and rapid collisional recombination 
and heating of the electrons d CHU > although we do see some signs of these 

effects, even in this regime. The ions typically equilibrate with Tj sa 1 K, which gives 
Ti « 3, so the ions are strongly coupled. 

The time scale for disorder-induced heating is the inverse of the ion plasma 
oscillation frequency, l/u) p i = y rriieo/nie 2 , which is on the order of 100 nanoseconds. 
Physically, this is the time for an ion to move about an interparticle spacing when 
accelerated by a typical Coulomb force of e 2 /Aireoa 2 . This time scale is also evident in 
molecular dynamics simulations of ion-ion thermalization |15lll6llT^ll20ll2~Tl 122112311^4] . 

For t delay S> 1/uipi, the ions have equilibrated and the thermal energy of the 
electrons begins to dominate the evolution of the plasma. Electrons contained in the 
potential created by the ions exert a pressure on the ions that causes an outward radial 

| As the number of electrons per Debye sphere (re -3 ) approaches unity, the Yukawa interaction ceases 
to accurately describe ion-ion interactions. For strontium plasmas studied here, this situation only 
occurs for the highest n e and lowest T e . It will be interesting to test Eq. |2]for these conditions. 
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acceleration. This was studied experimentally in jS] and theoretically by a variety of 
means in 17.. The experiments measured the final velocity that the ions acquired, 
which was approximately v ter minai ~ y/E e /rrii. With the imaging probe, we can now 
observe the expansion dynamics at much earlier times during the acceleration phase. 

As discussed in jS] and [TT|. a hydrodynamic model, which describes the plasma 
on length scales larger than Ad , shows that the pressure of the electron gas drives the 
expansion through an average force per ion of 

- -V(n e (r)fc B T e ) „rk B T e 
rii{r) af 

where the electron and ion density distributions are n e (r) w rij(r) = noiexp(— r 2 /2af). 
We assume thermal equilibrium for the electrons throughout the cloud |17|. 
The force leads to an average radial expansion velocity for the ions, 

V(r,t d elay) = ? ^—T^delay (4) 

rrnaf 

The velocity is correlated with position and increases linearly with time. This does 
not represent an increase in the random thermal velocity spread or temperature of 
the ions. Due to the large mass difference, thermalization of ions and electrons |15| is 
slow and occurs on a millisecond time scale. 

Equation^for the average ion velocity assumes a constant electron temperature. 
Actually, as the plasma expands, electrons will cool. This can be thought of in terms 
of energy conservation or adiabatic expansion. It is possible to describe the expansion 
with a Vlasov equation that includes the changing electron temperature. For an 
initial Gaussian density distribution, the equations can be solved analytically and 
the expansion preserves the Gaussian shape with a 1/s/e density radius given by 
°f (*) ~ °f (0) + [ k BT e {Q)/mi}t 2 delay [TTIEIIEEl. The ex P erimcnts involving absorption 
imaging of the plasma, however, have concentrated on the first few microseconds 
of the expansion when the plasma size and electron temperature have not changed 
significantly. Thus we can safely use Eq. 21 The effects of the expansion are evident 
in the radial velocity that manifests itself in Doppler broadening of the ion absorption 
spectrum. 

4. Doppler-Broadened Spectrum 

To obtain quantitative information from the plasma images, we relate the OD (Eq. 
to underlying physical parameters. Following Beer's law, the OD for a laser 
propagating along the z axis is 

OD(x,y)= J dzni(x,y,z)a[u,Ti(r)], (5) 

where rii(x, y, z) is the ion density, and a[v, Ti(r)] is the ion absorption cross section at 
the image beam frequency, v. The absorption cross section is a function of temperature 
due to Doppler broadening, and since we expect the temperature to vary with density, 
we allow a to vary with position. If we now integrate over x and y, or, in reality, sum 
over the image pixels multiplied by the pixel area, we get the spectrum 

S(u) = j dxdyOD(x,y) = J d 3 r m(r)a[v, T,(r)], (6) 
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Figure 4. Absorption spectra of 
ultracold neutral plasmas. We plot 
the integral of the optical depth (Eq. 
l6l . The frequency is with respect to 
a Doppler-frec absorption feature in a 
strontium discharge cell. Both spectra 
correspond to T e = 56 K and an initial 
peak plasma density of noi = 2 X 



K) 1 



Data are fit with Voigt 



profiles, and the increase in linewidth 
for longer t delay is clear. 



as a function of the image laser detuning§. As we vary the detuning, we obtain 
absorption spectra as shown in Fig. The rest of the paper will deal with 
the relationship between spectra such as these and the underlying temperature 
distributions of the ions. 

The absorption cross section for ions in a region described by a temperature T, 
is given by the Voigt profile 

{ , J 2 7e//l + 4(^)3 VWT) 6 



(7) 



where o\d(T) = ^/fcsT/mi/A is the Doppler width, and j e ff = 70 + Jiaser is 
the effective Lorentizian linewidth due to the natural linewidth of the transition, 
70 = 2-7T x 22 x 10 6 rad/s, and the laser linewidth, 72 Qser = 2n x (5 ± 2) x 10 6 
rad/s. The center frequency of the transition is vq = c/A, where A = 422 nm. The 
"three-star" symbol, 3* = 1, accounts for the equal distribution of ions in the doubly 
degenerate ground state and the linear polarization of the imaging light [221 ■ We fit 
the spectrum to a single Voigt profile with effective temperature T^ e ff, 

S{v)= J , d 3 rn i (r)a[i/,T i (r)]=JV i a(i/,T i , e// ), (8) 

where N% is the number of ions. We will explain below that T. e y/ is a good 
approximation of the average ion temperature in the plasma. 

5. Recent Results 

FigOJa) shows the evolution of T )6 j t for three different densities. The rapid increase in 
the temperature for t,i e iay < 300 ns is due to disorder-induced heating. This originates 

§ We can also fit OD(x, y) to a two dimensional Gaussian, as described in 1111 . and identify 
J dxdyOD(x,y) rj 2-KOi x criyOD max , where Oi x and <Ji y are the transverse sizes of the absorption 
profile, and OD max is the peak optical depth. This sometimes has signal-to-noise ratio advantages 
over integrating the entire image, but both approaches should give the same result. 
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Figure 5. The effective ion temperature, T i e ff, versus time after photoionization 
for initial electron temperature of T e = 2E e /3kB = 38 ± 6K and various plasma 
densities. (A) The data is plotted on absolute temperature and time scales. (B) 
The time is scaled by the inverse of the average plasma period, and e t f is scaled 
by T c . From 1501 . 



from conversion of Coluomb potential energy into kinetic energy as the ions evolve from 
a completely disordered state to one with some degree of spatial correlations. This 
was predicted in [23] , and observed in numerical simulations ^3 ESI EH and the first 
experimental studies with optical imaging 

The time scale of the heating is the inverse ion plasma frequency w^ 1 = 

y/niiEo/nie 2 w 100 ns. The energy scale is Tc = e 2 /47r£ a&B ~ 5K, where 
a = (47rrii/3) -1 / 3 is the Wigner-Seitz radius. Fig GJb) shows the same data with 
time scaled by the average inverse plasma frequency and temperature scaled by Tc- 
The three curves coincide very well in the time axis, but show slight deviation in 
the temperature axis. The deviation indicates the effects of electron screening of the 
ion-ion interaction. Electron screening is discussed in |30) . 

To quantitatively analyze the data in Figs. we must account for the effect of 
the inhomogeneous density distribution |31j . We expect the ion temperature to vary 
with density because global thermal equilibrium occurs on a hydrodynamic time scale, 
cr/v, which is on the order of tens of /lis, where v is the ion acoustic wave velocity. 
Local thermal equilibration occurs on a much faster time scale, ~ a; " 1 | 32 | . For a 
range of reasonable temperature distribution, numerical simulations |31| show that 
Ti,eff — (0.95 ± 0.05)Ti ave . Here Ti^ ave is the average ion temperature in the plasma. 
The extracted Ti iave and Eq. agree within our uncertainties of typically 0.2 K. 

If the effective ion temperature is measured for longer delay times, acceleration 
of the ions by the electron pressure also becomes important. The effect is small for 
the data presented here, but its effect can be modelled 123 EH]- The acceleration of 
the ions then becomes a sensitive probe of the electron temperature which evolves 
under the influence of recombination, continuum lowering, disorder-induced heating, 
and adiabatic expansion EH EE El EHj 

From the measured Ti^ ave , T e , and density distribution, we calculate the average 
Coulomb coupling constant for Debye-screend ions, 

Kavg = (eM-<r))e 2 /(4ire a(r)k B T have )}. (9) 

For all the data shown in Figs. Eland in pfll] . r£ is in the range of 1.7 to 2.5. With 
lower T e and higher density, r* atJe is slightly higher. The surprisingly small variation 
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Figure 6. Effective ion temperature obtained from different selected regions of 
the cloud for n 0i = (4 ± 2) X 10 9 cm" 3 and initial T e = 2E e /3k B = 38 ± 6K. In 
(b), we plot the visibility of the oscillation for the central probed region versus 
the averaged re for this region for all our data taken at different conditions. From 

1501 . 

in r* avg suggests that disorder- induced heating is a natural feedback mechanism that 
leads to equilibration just barely in the strongly coupled regime. 

Close inspection of Figs. |5] reveals that at the end of the disorder- induced heating 
phase, the ion temperature overshoots it equilibrium value before settling. This 
phenomenon is more evident in Fig. Eta), where Ti^ ave is calculated for an inner and 
outer region of the plasma image (p = \J x 1 + y 2 < 0.9a and p > 1.48tr respectively). 
Each selected annular region contains 1/3 of the ions and probes a region with 
significantly less variation in density than in the entire plasma. The region with 
lower density has lower Ti iave , supporting our hypothesis of local thermal equilibrium, 
but the oscillation is the most striking observation. 

Intuitively, one can explain this phenomenon as an oscillation of each ion in its 
local potential energy well. A simple calculation implies that the time for an ion to 
move an interparticle distance, when accelerated from rest by a force of e 2 / (Att£qo 2 ) 1 
varies as ~ w" 1 , and the observed oscillation occurs at 2u> pi . As expected from the 
density dependence of tu P i , the oscillation period is longer for the outer region where the 
average density is lower. This explains why averaging over the entire cloud obscures 
the oscillation; the motion dephases because of the variation in u P i. It is questionable 
whether the ion motion should be called an ion plasma oscillation or not, because 
there is probably no collective, or long range coherence to the motion. 

Kinetic energy oscillations at 2u> p i have been observed in molecular dynamics 
simulations of equilibrating strongly-coupled systems [20IE3E1- Calculations |3"3ll3"3) 
also show oscillations in the velocity autocorrelation function in equilibrium strongly- 
coupled one-component plasmas. To our knowledge, this is the first experimental 
observation of the phenomenon. 

Numerical results |2Q| suggest that the damping time for the oscillations is 
approximately ir/u> P i for Ti > 5, and that lower Ti leads to faster damping. Because 
our analysis still averages over the z-axis of the plasma, which introduces dephasing 
into the observed oscillations, it is difficult to comment on the damping of the 
oscillation. But we do find a correlation between the visibility of the oscillation and k, 
as shown in Fig.|Hfb), suggesting that electron screening plays a role in the damping. 
The visibility is defined as - Tf% f )/(TQ) + T^ ff ), where T%f f and T** f 
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are effective ion temperatures at the peak and dip of the oscillation. 

In conclusion, we have developed absorption imaging as a tool for studying 
ultracold neutral plasmas The absorption spectrum, after accounting for the 

inhomogeneous plasma density distribution, measures the average ion temperature. 
Equilibration of the plasma shows the effect of disorder-induced heating. The screening 
effect of electrons on the final ion temperature has been studied, and experiment and 
theory agree very well [23j- The ion kinetic energy also displays oscillations that are 
characteristic of equilibrating strongly coupled plasmas. 

Many future experiments suggest themselves. Some of the most interesting are 
investigating dynamics when the initial electron Coulomb coupling parameter is large 
and recombination and disorder-induced electron heating are expected to dominate 
the plasma evolution. Detailed examination of ion and electron thermalization at 
the border of the strongly coupled regime is also possible. Quantitative study of the 
kinetic energy oscillations is also necessary in order to develop a better understanding 
of aspects of the phenomenon such as damping. Improvements in imaging optics will 
also increase the image signal-to-noise ratio and allow study of features on the ion 
density distribution with ~ 10 /im experimental resolution. 

The fact that strontium ions scatter blue light so efficiently suggests one of the 
most exciting future directions for studies of ultracold nuetral plasmas. It may be 
possible to use laser light to cool and trap the ions in a strontium ultracold neutral 
plasma ^21; j us * as lasers are used to cool and trap neutral atoms. Confining the 
ions would also confine the electrons through Coulomb attraction. This would be an 
entirely new method of plasma confinement and may also lead to drastically lower 
plasma temperatures and stronger coupling for the ions. 

This research was supported by the Department of Energy Office of Fusion Energy 
Sciences, National Science Foundation, Office for Naval Research, Research Corpora- 
tion, Alfred P. Sloan Foundation, and David and Lucille Packard Foundation. 
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